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Abstract High angular resolution X-ray imaging is always demanded by 
astrophysics and solar physics, which can be realized by coded-mask imag- 
ing with very long mask-detector distance in principle. Previously the 
diffraction-interference effect has been thought to degrade coded-mask imag- 
ing performance dramatically at low energy end with very long mask- 
detector distance. In this work the diffraction-interference effect is described 
with numerical calculations, and the diffraction-interference cross correla- 
tion reconstruction method (DICC) is developed in order to overcome the 
imaging performance degradation. Based on the DICC, a super-high angular 
resolution principle (SHARP) for coded-mask X-ray imaging is proposed. 
The feasibility of coded mask imaging beyond the diffraction limit of single 
pinhole is demonstrated with simulations. With the specification that the 
mask element size of 50 x 50 /im 2 and the mask-detector distance of 50 m, 
the achieved angular resolution is 0.32 arcsec above about 10 keV, and 0.36 
arcsec at 1.24 keV (A = 1 nm) where diffraction can not be neglected. The 
on-axis source location accuracy is better than 0.02 arcsec. Potential appli- 
cations for solar observations and wide-field X-ray monitors are also shortly 
discussed. 
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1 INTRODUCTION 

High angular resolution X-ray imaging is always demanded by astrophysics and solar 
physics, for example to study black holes near event horizon, relativistic jets of super mas- 
sive black holes, as well as solar flares and coronal activities. So far the best imaging tech- 
nology for X-ray observation is realized by grazing incidence reflection fAschenbach ll985f , 
which provides a very good angular resolution, for example down to 0.5 arcsec as in the 
case of the Chandra X-ray Observatory (Weisskopf et al. 2000]). The diffraction limit of 
Chandra is about 8 miliarcsec at 6 keV; however slope errors and other surface irregu- 
larities that affect the angle of reflection prevent the angular resolution better than 0.5 
arcsec. Further more, the grazing incidence reflection is limited by the working energy 
band, which hardly exceeds tens of keV. A possible way to improve the angular resolution 
beyond Chandra's requires a different technology, for example the diffractive-refractive 
X-ray optics proposed by Gorenstein (j2007|) or the Fourier- Transform imaging by Prince 
et al. (1988]). As an alternative, a super-high angular resolution principle (SHARP) for 
a coded-mask X-ray imaging telescope concept is proposed here. 

The coded mask imaging technology, which has been reviewed by Zand (|1996p . is 
widely applied for X-ray observations, for example the INTEGRAL mission (Winkler 
et al. 1211051) of ESA and the SWIFT mission (Gehrels I21TO)) of NASA. Both missions 
have angular resolution in tens of arcmin, since the distances between the masks and the 
detectors are limited by the dimensions of the satellites. In principle super-high angular 
resolution better than Chandra can be achieved by a coded-mask telescope with sub- 
millimeter size of mask pinholes and long mask-detector distances, which can be tens, 
even hundreds of meters with the formation-flying technology (for example proposed for 
XEUS (ESA l200T|) of ESA) or the mast technology (for example applied in Polar mission 
of NASA). The X-ray diffraction effect in a coded-mask system has been commonly 
thought to be negligible due to high photon energies in hard X-ray and gamma-ray range 
as well as relatively small mask-detector distance. Whereas for coded-mask telescope 
with a super high angular resolution in the milli-arc-second range, the X-ray diffraction 
effect could not be neglected, especially at low energy end (in the range of several keV), 
which might degrade the imaging performance remarkably, as pointed out by Prince et 
al. (|1988p and Skinner (j2004|) . In this work, we study the diffraction effect in SHARP 
with numerical computations, and propose a diffraction- interference cross correlation 
reconstruction method and demonstrate the feasibility of coded mask imaging beyond 
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the diffraction limit of single pinhole. A potential application for solar observations is 
also shortly discussed. 

2 THE CODED-MASK PRINCIPLE AND 

DIFFRACTION-INTERFERENCE RECONSTRUCTION 

2.1 The coded- mask principle and SHARP 

In the limit of geometrical optics, i.e., the diffraction effect is not considered, the ba- 
sic concept of coded-mask imaging is shown in Fig. [T] The coded-mask camera has a 
mask on top of a position-sensitive detector with a mask-detector distance D. Two point 
sources project the mask pattern onto the detector plane. The shift and the strength of 
projections encode the position and the flux of sources separately. The detection of the 
X-ray flux can be described with Equ. [T] (Fenimore et al. I1978[) , 

P = O * M + N noise , (1) 

where O is the X-ray flux spatial distribution, M is the encoding pattern of mask and 
* means cross correlation. For the commonly applied cross correlation method, a matrix 
G is used for reconstruction, as shown in Equ. [2l 

O 1 = P*G = 0*(M*G) + N noisc * G, (2) 

where O' is the estimation of X-ray flux spatial distribution, if M*G is a 6 function. The 
encoding pattern M is often chosen so that its auto-correlation function is a 8 function. 
Therefore G = M is commonly applied. The width of the Point Spread Function (PSF) 
obtained with Equ. [His defined as A, = ^g- (d m is the size of the mask element). The 
cyclic optimal coded-mask configuration ( Zand [T996P is also employed here. The coded 
mask or the mask pattern consisted of four basic patterns is about four times the size 
of the detector, and one basic pattern cyclically extends about twice along the width 
and length directions of the mask. In the next sections, simulations are done with the 
following assumptions: the basic pattern is totally random; the detector pixel size is equal 
to the mask element size, both of which are in square shape; the mask and the detector 
are assumed to be idealized, i.e., the mask has no thickness, the opaque mask elements 
block the X-ray completely and one incident event can only be recorded in one detector 
pixel. 

In this work, d m — 50 /im and D = 50 m are chosen to be the basic parameters 
for the simulations of SHARP, which means Aj = 0.2 arcsec in sub-arcsec range. Since 
the simulations employ two dimensional FFT (please refer to Sec. I2.2p which requires 
tremendous computing resources, the basic pattern simulated contains only 200 x 200 
elements, i.e., the reconstructed image contains the same pixels and each pixel is 0.2 x 
0.2 arcsec 2 . The fully coded Field-Of-View (FOV) is 40 x 40 arcsec 2 . 
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Fig. 1 The basic concept of coded-mask imaging, in the limit of geometrical 
optics, i.e., the diffraction effect is negligible. Two point sources illuminate a 
position-sensitive detector plane through a mask made of many square pinholes. 
The detector plane thus records two projections of the mask pattern. The shift 
and the strength of projections encode the position and the flux of the sources 
separately. The PSF width is defined as Aj = 



2.2 The diffraction-interference cross correlation method 

As shown in Fig. [21 A is the mask modulation function, which is 1 or with idealized mask, 
r is a vector from point (xo, j/o) on the mask to the (x, y) on the detector, k is the wave 
vector. Since there are many pinholes in the mask, both diffraction and interference will 
take place on the detector plane. The diffraction-interference is described with Fresnel- 
Huygens principle (Lindsev ll978p as, 

poo poo ikr 

Eii(x,y) = C / / Ei(x ,yo)A(x a ,y ) dx dy , (3) 

J oo J oo ^ 

where Ei is the amplitude distribution of incident photons, which is described with the 
plane wave in our case as Ei — A4 * e 1 . The flux distribution is proportional to Ef. 
C = —i/X is constant. In our case the dimension of the mask is far less than the mask- 
detector distance D, we have r ~ D + - +x ° +v +Vo^ 2xx o-'2vvo ^ Therefore Equ. [3] becomes 

/>oo />oo 

E R (x,y)=C 2 / Me- 2i *V* Xo+ fyy°Ux dy , (4) 

J DO J CO 

where C 2 = f^e ikD 'e^( x '' f x = x/XD, f y = y/XD, M = 
Ei(xo, yo)A(xo, yo)e2v( x ° +Vo ' > ■ Equ.[4]is the Fourier transform approximation of Fresnel- 
Huygens principle, which can be simulated with an FFT algorithm. The conditions for 
the validity of this approximation are discussed by Lindsey (|1978p . and are satisfied in 
our simulations. 

The X-ray diffraction can not be neglected at the low energy end as in Fig. [3l which 
shows the simulated flux distribution on the detector plane with a 50 /im square pinhole 
in front of it (50 m away) for 1.24 keV photons (A — 1 nm). The flux spreads much larger 
than 50 /urn (each pixel in Fig. [3] is 50 x 50 /im 2 ) due to the diffraction effect. 
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Fig. 2 Monochromatic X-ray with wave vector k illuminating onto a mask in 
the XQ-yo plane produces a diffraction pattern in the x-y plane (Lindsey I1978|) . 



If the reconstruction matrix G in the cross correlation method is still chosen to be 
equal to M (so-called geometrical reconstruction matrix, which is indicated as Go) at the 
low energy end, the diffraction effect would degrade the reconstructed image remarkably. 
Fig. [4£i shows the reconstructed image of a monochromatic on-axis point source (A = 
1 mm) obtained by geometrical matrix. The reconstructed image is normalized so that 
the total integral flux is one photon, and all the reconstructed images in this work are 
normalized in the same way if not specified. The point source is hardly identified with 
the Fresnel diffraction stripes clearly observed. Since the X-ray probability wave through 
different mask open elements interferes with each other, the mask spatial information 
should be encoded into the diffraction-interference pattern on the detector. We then 
choose the reconstruction matrix G at certain energy as the normal incident diffraction- 
interference pattern (so-called diffraction- interference matrix, which is indicated as Gdi) 
on the detector plane. Fig. UJd shows the reconstructed image obtained with Gdi of 
the same source in Fig. [4^,, and the point source can be clearly identified although slight 
Fresnel stripes still exist. This new reconstruction method is called diffraction- interference 
cross correlation method (DICC), which can be described with Equ.[5j 



Of = P * Gdi = O * (M * Gdi) + N noisc * Gdi- 



(5) 



Gdi can be obtained by numerical simulation as shown in Equ. 2] or from actual mea- 
surements. The angular resolution of DICC is beyond the single pinhole (the single mask 
open element) diffraction limit as indicated in Fig. 21 
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Fig. 3 The simulated flux distribution on the detector plane for 1.24 keV pho- 
tons (A = 1 nm) with a 50 fim square pinhole in front of it (50 m away). The flux 
spreads much larger than 50 /im caused by diffraction as one pixel is 50 x 50 /im 2 . 
Therefore this diffraction effect must be included when reconstructing the image 
of the incoming X-ray with cross-correlation effect. 
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Fig. 4 The reconstructed images of a monochromatic (A = 1 nm) on-axis 
point source: (a) the cross-correlation with the Go, the diffraction degrades 
the reconstruction dramatically, because the diffraction pattern shown in Fig. 3 
has not been taken into consideration; (b) with DICC, the point source can be 
clearly identified although slight Fresnel diffraction stripes still exist. The DICC 
can thus largely overcome the degradation of the imaging performance due to 
the diffraction of each pinhole. 



3 THE SIMULATION RESULTS 

We have simulated the reconstructed images at several energies from 1.24 keV (A = 1 
nm) to 6.2 keV (A = 0.2 nm), where DICC is valid. To illustrate the main properties of 
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this method, we mainly discuss the imaging performance of DICC around 1 nm, which 
is also compared with that in the limit of geometrical optics. 

3.1 Angular resolution 

We calculate the reconstructed image of two point sources with the same flux (Poisson 
fluctuations are not considered here). The angular resolution is defined as the minimum 
angular distance between these two sources, which can be separated by the 50% contour 
line of the maximum flux on the reconstructed image. Fig. [5^ shows the reconstructed 
images in the limit of geometrical optics for two point sources with angular distance 
1.5Ai = 0.3 arcsec; the two sources can not be separated. However in FigEb two sources 
with angular distance 1.6Aj = 0.32 arcsec can just be separated. Therefore the angular 
resolution is 0.32 arcsec in the limit of geometrical optics, i.e., at high energy end (above 
about 10 keV) with diffraction negligible. With the same procedure, the angular resolu- 
tion of DICC for A = 1 nm photons is obtained as 0.36 arcsec, i.e., the two sources can 
not be separated in Fig. [St, but can be separated in Fig. [SJi. This means the diffraction 
effect degrades the angular resolution by only about 12.5%, a significant improved com- 
pared to that in Fig. 3. We also calculated the reconstructed images with DICC for A = 1 
nm point sources at different locations inside the FOV, but no observable reconstruction 
degradation is found compared with the on-axis source case. 

3.2 Source location accuracy 

The reconstructed images for an on-axis point source is simulated 10000 times with 
1000 photons recorded (Poisson fluctuations are considered). The source location of each 
simulation is calculated by fitting the reconstructed image to the PSF. The calculated 
source location distribution with DICC for A = 1 nm is shown in Fig. [5^, which is 
concentrated within O.lA^ = 0.02 arcsec. The same simulation is also done in the limit 
of geometrical optics, which is shown in Fig. [6Jd. Therefore the source location accuracy 
is better than 0.1 Aj = 0.02 arcsec for 1000 detected photons. 

3.3 System constraints 

Since the diffraction pattern on the detector plane is energy dependent, the diffraction 
matrix G is also energy dependent. In practice the detector can not identify two photons 
with an energy difference within the detector's spectroscopy resolution, i.e., in practice 
photons within a small energy band must share a common diffraction matrix. Fig.[7]shows 
the reconstructed images of several on-axis point sources with wavelengths around 1 nm, 
which all apply the diffraction matrix of 1 nm photons (marked as G\ ) . The point source 
can be clearly identified. Obviously the imaging qualities are degraded by applying the 
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Fig. 5 The angular resolution of the simulated system: (a) and (b) are in the 
limit of geometrical optics; (c) and (d) are obtained with DICC for A = 1 nm 
sources. In (a) and (c) two sources can not be separated. However in (b) and 
(d) two sources can be separated. The angular resolution is 0.32 arcsec in the 
limit of geometrical optics and degraded to 0.36 arcsec for A = 1 nm sources. 



Table 1 The Degradation of the Angular Resolution with Shared 
Reconstruction Matrix 



Wavelength (nm) 0.95 0.97 0.99 1 1.01 1.03 1.05 

Angular resolution (arcsec) 0.39 0.38 0.37 0.36 0.37 0.39 0.41 



G\ instead of their own diffraction matrixes. Tab.[T]gives the angular resolution at several 
wavelengths around 1 nm by applying G\ as reconstruction matrix instead of their own. 
For the 1.05 nm on-axis source, the angular resolution is degraded to 2.05Ai = 0.41 arcsec 
by applying d. 

Based on the detector energy resolution, the system constraint is considered. 
Assuming the detector pixel size dd, the difference Sd of the main diffraction stripe 
diameters of photons with wavelength difference <5A is 



_ 5\*D 

dm 



(6) 
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Fig. 6 The calculated source location distribution of the reconstructed images 
for an on-axis point source, which is simulated 10000 times with 1000 photons 
recorded: (a) with DICC for A = 1 nm photons; (b) in the limit of geometrical 
optics. The source location accuracy is better than O.IA^ = 0.02 arcsec. 



If the spectroscopy resolution of the detector is 5X, spatially the detector does not 
need to distinguish Sd for photons with wavelength difference SX, 

Sd < f * d d . (7) 

Here / > 1 is a constant and depends upon the mask open fraction. For a 50% open 
mask, our simulations indicates that 1 < / < 2. For a mask fraction much less than 50%, 
/ > 2 is possible. Then with known basic parameters, i.e., d m and D, the requirement 
for detector spectroscopy resolution is 

SX < f * A, * d d . (8) 

In practice d rn and D are chosen to match the detector performance, which means 

dm > ~~F~~/C~) (9) 

J * Ai 



^ = ^<^< 7XW . (10) 



or 

dm < dd_ SX 
A, " A, " (A,) 2 /' 

Therefore the larger the / is, the more compact the system can be. Here / = 1 is 
chosen for convenience, with which SX = 0.05 nm, d m — 50 fjxa and D = 50 m satisfies 
the system constraint. The spectroscopy resolution 0.05 nm requires about 120 eV energy 
resolution @A = 1 nm, which can be provided by modern silicon imaging detectors like 
CCD or DEPFET (Struder [20001). 

4 DISCUSSION AND CONCLUSION 

By applying DICC, the system angular resolution is no more limited by single pinhole 
diffraction limit at low energy end, but limited by detector spectroscopy performance. 
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Fig. 7 The reconstructed images of (a) 0.9 nm (b) 0.95 nm (c) 1.05 nm and 
(d) 1.1 nm photons, which all apply the same diffraction- interference matrix G\ 
(the diffraction-interference matrix of 1 nm photons). 



Tab. [2] shows a brief comparison of SHARP with several X-ray astronomy missions. The 
advantage of SHARP is quite obvious: excellent angular resolution and easily fabricated 
mask which does not require high accuracy mechanical fabrication of X-ray grazing lens. 
However the detector effective area of a coded-mask telescope is smaller than the its actual 
area due to the encoding pattern. Further more each detector pixel receives background 
photons from all over the FOV. Therefore the sensitivity is limited by its large background 
and small effective area compared with that in the grazing incidence reflection case. 

Current detector technology can provide even better energy resolution than 0.05 nm 
@1 nm, which means the angular resolution discussed above can be improved further 
by increasing the mask-detector distance D as indicated in Equ. [TU] or by decreasing 
the mask element size as indicated in Equ.[9l For example for a detector with 0.002 nm 
spectroscopy resolution (5 eV @ 1 keV) , d m = 40 /im and D = 800 m satisfies the system 
constraint with A, = 0.01 arcsec. 

A potential application of SHARP will be the solar observation at 1-100 keV en- 
ergy band with sub-arcsec angular resolution. The main objectives may be the coronal 
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SHARP Chandra Integra^ RHESSI Hinode 

(Solar B 



Mission 



Imaging 
technol- 
ogy 

Energy 
band 
Energy 
resolution 



Angular 
resolution 



Coded- Focusing Coded- Rotation Focusing 

mask mask Modulation 

& Fourier 
Transform 

1-100 keV Up to 10 15 keV-10 3-150 

MeV ke\|] 
9 % @ 100 < 1 kevQ 
eV@5.9 about 148 keV 
ke\|] eV@5.9 
ke\|] 

0.32-0.36 0.5 arcseqj 12 arcmin 2.26 arc- 2 arcsec 



133 



keV 
ACIS: 



About to 6 
keV 

< 280 eV 
Q5.9 ke\i 



0.523ke 



a IBIS onboard 

b XRT onboard 

c The upper layer 

d Take DEPFET as focal plane detector (Treis et al. |20"06)> 

e CCD S3 on board 

s Below 100 keV 

9 Calculated based on the ENC < 30 el. 

h On-axis point source 

' The collimator # 1 

3 68% source flux within 2 arcsec 



mass ejection (CME) and solar flares, including fine structures and evolutions of the 
solar flares, nonlinear solar flare dynamics, solar particle acceleration mechanism et al. 
Therefore applications of SHARP are foreseen to make significant progress on the study 
of solar high energy explosive events and the space weather forecast model. Another po- 
tential application is to make wide-field X-ray monitors with SHARP; its sub- arc angular 
resolution may allow the counterparts of gamma-ray bursts, X-ray bursts, black hole and 
neutron star transients to be identified without the requirement of subsequent follow-up 
observations of focusing X-ray telescopes. 
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